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Abstract 

Kinetic-transport equations are, by now, standard models to describe the dynamics of popula¬ 
tions of bacteria moving by run-and-tumble. Experimental observations show that bacteria increase 
their run duration when encountering an increasing gradient of chemotactic molecules. This led to 
a first class of models which heuristically include tumbling frequencies depending on the path-wise 
gradient of chemotactic signal. 

More recently, the biochemical pathways regulating the flagellar motors were uncovered. This 
knowledge gave rise to a second class of kinetic-transport equations, that takes into account an 
intra-cellular molecular content and which relates the tumbling frequency to this information. It 
turns out that the tumbling frequency depends on the chemotactic signal, and not on its gradient. 

For these two classes of models, macroscopic equations of Keller-Segel type, have been derived 
using diffusion or hyperbolic rescaling. We complete this program by showing how the first class of 
equations can be derived from the second class with molecular content after appropriate rescaling. 
The main difficulty is to explain why the path-wise gradient of chemotactic signal can arise in this 
asymptotic process. 

Randomness of receptor methylation events can be included, and our approach can be used to 
compute the tumbling frequency in presence of such a noise. 


Key words: kinetic-transport equations; chemotaxis; asymptotic analysis; run and tumble; biochem¬ 
ical pathway; 
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1 Introduction 

Two classes of kinetic-transport equations have been proposed to describe, at the cell scale, the 
movement of bacteria by ‘run and tumble’ in a given external effective signal M{x, t), usually related 
to the extra-cellular chemo-attractant concentration 5 by a relation of the type M = mo + ln(5). 
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The simplest class is for the probability p{x,v,t) to find a bacteria at location x G and with 
velocity v ^ V {a. smooth bounded subset of one can choose the unit ball to fix idea ). The 
evolution of this probability is given by a Boltzmann type equation 

dtp + V-VxP = T[DtM]{p), X gR'’-, V £ V, t >0, (1) 

with the path-wise gradient of M defined as 

DfM = dtM + V ■ VxM ( 2 ) 

and T a tumbling kernel which typically takes the form 

T[DtM]{p) = [ [T{DtM{x,v',t),v,v')p{x,v',t)-T{DtM{x,v,t),v',v)p{x,v,t)] dv'. (3) 

Jv 

Such equations, with T depending on M or DtM, were used intensively to model bacterial chemotaxis, 
possibly with M connected to the cell density, as a result of chemoattractant release by bacteria. They 
were first introduced in [T3] and the Keller-Segel drift-diffusion system was subsequently derived |15j 
in the diffusion limit; surprisingly, with a kernel T depending on M and not on its gradient, and in 
opposition to the Keller-Segel system which solutions blow-up for large mass, it was proved that the 
solutions exist globally mm- However, experiments show that bacteria as E.coli extend their runs 
when feeling an increasing concentration of chemoattractant and this led to study tumbling kernels T 
that depend on DtM, see mm- The nonlinear theory is then more difficult (see [3] and the references 
therein) and blow-up can occur in finite time [2]. These models with T depending on DtM are able 
to explain the experimental observation of traveling pulses of bacteria, which cannot be done when 
T only depends on M itself, see [niiiT]. Also, departing from this kinetic-transport equation, it is 
possible to rescale it and study the diffusion and hyperbolic limit as in [3 muni EH Em i- When T 
undergoes stiff dependency on DtM, the hyperbolic limit is singular and the analysis is particularly 
delicate mi- 


More elaborated kinetic models have been proposed recently that incorporate intracellular chemo- 
sensory system. In the simplest description of the biochemical pathways, they use a single additional 
variable m > 0, which represents the intracellular methylation level. Then, the kinetic-transport 
equation is written for the probability density function p{x,v,m,t) of bacteria at time t, position 
X gR'^, moving at velocity v G V and methylation level m > 0 


dtp + v- VxP + dm[f{rn,M)p] = Q[m,M]{p), 
p{x, v,m = 0,t) = 0. 


( 4 ) 


The intracellular adaptation dynamics is described by the reaction rate /(•) for which we assume 
f{m = 0,M) > 0, which allows us to pose the boundary condition at m = 0. The tumbling term 
Q[m, M]{p) is 


Q[m,M]{p)= / [X{m,M,v,v')p{t,x,v',m) — X{m,M,v',v)p{t,x,v,m)] dv', (5) 

Jv 

where X{m,M,v,v') denotes the methylation dependent tumbling frequency from v' to v, in other 
words the response of the cell depending on its environment and internal state. We borrow this for¬ 
malism from da EH] even though this type of models, involving more general signal transduction, can 
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be traced back to laBEi 121) . The authors in [3 Cl 0 [H EB [22] developed the asymptotic theory 
which allows to recover, in the diffusion and in the hyperbolic limits, macroscopic equations where the 
variables are only {x,t) as the Keller-Segel system, or {x,m,t) for structured Keller-Segel models. 

In the program of establishing the relations between these pieces of the model hierarchy for bacterial 
population motion, a derivation is missing: how are related these two classes of kinetic models 0-0 
and 0-0? 

Our goal is to show how, assuming fast adaptation and stiff response, the methylation level is at 
equilibrium with the external signal represented by M, and the equation Q can be derived from Q. 
In particular we aim at computing the bulk tumbling kernel T{DtM,v,v') from the methylation 
dependent kernel X{m, M,v,v'), a statement we give in the next section. Two difficulties arise here: 
one is to infer the proper rescaling in the kinetic equations, the second is to carry-out the mathematical 
analysis for singular limits. Our approach allows us to also include noise resulting from random 
receptor-methylation and demethylation events. The proof of the formula for T is given in sections 
and we show that a direct use of the variable m is not enough to produce the formula and that 
a new variable is needed, which zooms on the intra- and extra-cellular methylation equilibrium. We 
conclude by relating our notations to a more physically based description of the same model where 
the cell receptors activity is used in the model parameters, see section 

To keep simplicity, we assume that the external signal function M(x, t) is given and smooth. There¬ 
fore questions of existence and blow-up are not considered here. 


2 Fast adaptation, stiff response 

Assumptions. For our mathematical derivation, we introduce a small parameter e which acts both 
as a fast time scale for external signal transduction and as a stiffness parameter for the response in 
terms of tumbling rate. We assume moreover that the reaction rate / only depends on the difference 
m — M, in accordance with the physical models that we recall in Section Therefore, we rescale 
equation 0-0 as 


dtPe + v ■ VxPe + ldm(^f{m - M)p^ = Q,[m,M](p^), 

t;, m = 0, t) = 0, 

with the tumbling kernel 


Qe[m,M]{p^) = 


Iv 


A 


m — M 


/ N A Z™- “ \ / N 

,v ,m,t) — - ,v , vjp^[x, V, m, t) 


dv'. 


We complete this equation with an initial data > 0 which satisfies 


(1 -|- w?)p™'^{x, V, m) dx dv dm < oo. 


*xVxI 


:= / G x V). 


( 6 ) 

(7) 

( 8 ) 

(9) 
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Also, we are going to use several assumptions for the functions M, f and A. We assume they are 
as smooth as necessary and that for some constants m±, g±, A±, 

0 < m_ < M{x,t) < m+, M G x [0, oo)), (10) 

f{y) = -yGiy), with G G C'b(M), 0 <g- < G{y) < g+, (11) 

dyA(^y,v,v'^ <0, 0 < A_ < A(^y,v,v'^ < A+. (12) 

Various scalings have been proposed for kinetic equation and the closest, but still different seems to 
be the high field limit [1]. In m, still other scalings or limits are studied. 


The main result. With these assumptions, we are going to show that as e vanishes, we recover the 
simpler model Q-Q as a limit of ([^. 


Theorem 2.1 (Derivation of the kinetic equation) We make the assumptions (U)-® on the 
initial data, and (|10[)^( 12) on the coefficients. Let p^ he the solution to (§. Then, for all T > 0, 
is bounded in L°°(|o(T] x x V) and 


Pe := / Pedm po in L“([0,T] x x V)-weak-k 
Jr ^^0 

and po satisfies equation with 

Furthermore, we have po := f^podm with po the weak limit (in measures, see (dH); of pe which is 
given by 

Po{x, V, m, t) = po{x, V, t)5{m = M{x, t)). 


Before we prove this theorem in the next sections, we present a variant of this result. 

Internal noise. Due to random receptor-methylation and demethylation events, some internal noise 
can be observed in E. coli chemotaxis and we can model it by adding a diffusion term in m [6]. The 
model is as follows: 

dtPe + v ■Vxp,, + ^dm(^f{m-M)p,,^ = + Q,[m,M](pe), (13) 

with the no-flux boundary condition that now reads 

f{—M)pfix,v,m,t) — e^dmPe{x,v,m,t) = d, at m = 0. (14) 


Theorem 2.2 (Limit with noise) With the assumptions and notations of Theorem 2.1, the same 


conclusions hold for the solution p^ of (13), with the same expression for pQ and 


T[u,v,v') = [ A(y,v,v')e ^ 

V 27r 
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A priori bounds and principle of the proof. Before we explain the derivation of the formula 
stated in these theorems, let us make some observations which explain the difficulty. Because we 
assume that M(x,t) is given, we handle a linear equation for which existence and uniqueness of weak 
solutions is well established. The nonlinear case, when the chemoattractant concentration giving rise 
to M is coupled to p^, can also be treated, see m- In particular we will make use of the uniform 
estimates (see Section]^ 


'xl/xl 


Pe{x, V, m, t) dx du dm = 



'xFxI 


p^'^\x,v,m) dxdv dm, Vt > 0, (15) 


Pe{x,v,t) < \\p"’^\x,v)\\ooe^\ Vt > 0, (16) 

where C is a nonnegative constant. From these bounds, we conclude that we can extract subse¬ 
quences (but to simplify the notations we ignore this subsequence) which converge as mentioned in 
the theorems. 

Passing to the limit in the equation on (with or without noise) gives us 


dm(^f{rn - M)po) = 0. 


This tells us that f {m — M^po = 0 (it is constant and po is integrable). Because, with assumption 
f{m — M) vanishes only for m — M = 0, we conclude that po is a Dirac mass at m = M, hence the 


expression of po in Theorems 2.1 and |2.2[ 

However this information is not enough to pass to limit in the equation on p^ obtained integrating 
in m equation (1^ or (13), that is 


dtPe + V-Vj,Pe 


Qe[m, M](pe) dm. 


Indeed, in the right hand side, the product v'jp^{x, v', m, t) is, in the limit, a discontinuity 

multiplied by a Dirac mass. For this reason, we have to rescale in m in order to evaluate this limit, 
which we do in the next section. 


3 The change of variable 


To get a more accurate view of the convergence of p^ to a Dirac mass in m, and following [7], we 
introduce a blow-up variable around m = M. We set 


y = 


m — M 
e 


q^{x,v,y,t) = epe{x,v,m,t) 


(17) 


so that 


qe{x,v,t) := / qe{x,v,y,t)dy = / pe{x,v,m,t)dm = pe{x,v,t). 


(18) 


Because of these identities, our statements will equivalently be on q^ and will go through 
of g'e rather than p^ itself. 

Also notice that the bounds in (15), (16) also hold true for q^ and q^ and allow us to take 


the analysis 
weak limits. 
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(i) Without noise. The equation for q^[t,x,v,y) is written, using the definition in Q, 


dtqe + V ■ V^qe -^DtMdyq^ + ^dy{f{ey)qe) 


Iv 


^(y,'i^,v''^Qe{x,v',y,t) - A(y,v',v^q^{x,v,y,t) 


du^ 


From we can write f{ey) = eyG{ey) and the above equation becomes 
dtqe + V ■ VxQe -^DtMdyqe - ^dy{yG{ey)qe) 


Iv 


^(y,x,v''^q^{x,v',y,t) - A(y,v',v^q^{x,v,y,t) 


(19) 


du^ 


Because q^ is a bounded measure on x F x M'*' x (0, T), for all T > 0, as e —>■ 0, has a weak limit 
go in the sense of measure (again after extraction) and the above equation gives, in the distributional 
sense, 

dy{yG{0)qo + DtM{S)qo) 

From this, we infer that 

qo{t,x,v,y) = qo{t,x,v)6(y = - 
This information is useful provided we can establish that 


This step is postponed to Section and involves a control of the tail for large values of m. 
We may also integrate equation (|19l) with respect to y and find in the limit 


= 0. 

(20) 

DtM{S)\ 

G(0) )■ 

(21) 

c- lim qeix, v, t). 
e-S>0 

(22) 


dtqo + V ■ V^qo = 


Iv 


A^_m 

G(0) ^ 


■i;,u')godu'-A(^- 


DtM , 


(23) 


where D[M{S) is the total derivative, as in Q, but in the direction v' and where go represents 
qo{x,v',t). Finally, for any smooth test function (f), we have from the change of variable m ^ y = 
(m — M)/e, 

/ p,(a;,u,m,t)(/>(m)dm = / ge(a;, u, y, t)(/>(M + ey) dygo(®, u, t)(/)(M), 


where we use (21). This gives the limiting expression of po in Theorem 2.1 


These are the results stated in Theorem 2.1, if we can establish the relation go = po as stated in (22), 
which we do later. 


(ii) With internal noise. Similarly, after introducing the new variables as in (0, the equation ( [l^ 
for q^{t,x,v,y) writes 

dtqe + V ■ VxQe --^DtMdyqe + ^5y(/(ey)g,) 


= + 


Iv 


^(y,v,v''^qe{x,v',y,t) - A(y,v',v^qe{x,v,y,t) 


dv'. 
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From this equation becomes 


dtqe + V ■ VxQe -^DtMdyqe - ^dy{yG{ey)qe) 


= l^lyqe + 


'v 


^{y,v,v''^q^{x,v',y,t) - A(^y,v',v^q^{x,v,y,t) 


du'. 


In the limit e —)• 0, the above equation converges to, in the sense of distributions, 

dy{yG{0)qo + DtMqo) = -dyyqo, 

which shows that 

qoix,v,y,t) = qo{x,v,t) 2 {y am) ^ 


a useful information, still assuming we have proved the relation (22) for go- 


We conclude as before. After integration of (24) with respect to y, passing to the limit e 
find 

dtqo + v ■Vj,qQ = j - T^DtM^v.v'^qo udv', 

T(DtM,v,v''\=J~^^l K{y,v,v')e~~^'^^~^^ dy. 

\ ^ ^ Jm. 


The Theorem 2.2 is also proved, n 


(24) 


(25) 
0, we 

(26) 


4 A priori bounds 

We now establish the various estimates which justify that we can pass to the limit as indicated in 
Section and thus we prove the 


Lemma 4.1 We make the assumptions of Theorem \2.l\ then the condition (22) holds and for some 
constant which depends on ff y^g'“dydxdu and || Af||^i,oo(]gdxR+); have 


y‘^q^{t) dydxdv < G{q^^\M). 


Consequently, converges weakly in the sense of measure towards go and 


(i) for q^ a solution to (19), go is given by (21) with go weak solution of (23), 


(a) for g^ a solution to (24), then, go is given by (25) with go weak solution of (26) 


Proof. We only consider the case (i) without noise, the case (ii) is obtained by the same token. We 
first prove some estimates which imply weak convergence. Then, we pass to the limit in the equation 
satisfied by g^. 


bound. For completness, we recall that equation (19) is positivity preserving and conservative. It 


follows the uniform, in e, bound for g^ in L^, see (15). 


L°° bound on g^. We use the notation (18) for g^. Arguing in the spirit of uni 13120 ]), we first prove 
the uniform L°° bound on g^. 
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Integrating (19) with respect to y, from the bound (12) and the nonnegativity of q^, we get 


dtQe + V ■ j q^ du'. 


Then, using the method of characteristics, we have dtqe{t, x + vt) < A+ f qe{x + vt, v', t) du', which 
implies after integration 

qe{t,x,v) — q^'^'^{x — vt) < \+ / / q^ix — vs,v',t — s) dv'ds. 


Taking the supremum in x, v, we find 


+ A- 


1^1 f\ms)\\c 

Jo 


,ds, 


and using Gronwall’s inequality, we find the estimate in (16). 


Control on the tail in m. In order to prove the condition (22), we need to ensure that there is no 
mass loss at infinity in m. To do so, we multiply both sides of (19) by and integrate by parts with 
respect to x, v, and y. This yields 

— JJ y"^ q^ dy dx dv -\— JJ y‘^G{ey)qedy dx dv -\— JJ yDtMq^dy dx dv = 0. 

Using the Cauchy-Schwarz inequality, we deduce 
d 


dt 


JJ y‘^qe dydxduH— JJ y‘^G{€y)qf:dy dx dv <- JJ y^G{ey)qedy dx dv 


1 

+ - 
e 


{DtM)"^ f dy dx du. 


By assumption (10), DtM is bounded in L°°([0,r] x x V). From assumption © and the mass 
conservation, the last integral of the right hand side is uniformly bounded by a constant denoted by 
C > 0. Then from assumption (|11[), we have. 


d 


9 - 


^ I I yqedydxdv + — 


2 , , , C 

y qedydxdv< 

e 


From the Gronwall Lemma, we deduce the bound for all f > 0, 

y‘^qe{t) dy da; du < JJ y‘^q™' dy dx dv + 

which implies a uniform bound on ff y‘^q^{t) dydxdu. 
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Passing to the limit. From the bound above, we deduce that, we can extract a subsequence which 
converges weakly in measure q^ q^ and such that q^ % in L°°-weak-k. Then we can pass to the 
limit in the sense of distribution in equation (19) and deduce that the limit qo satisfies equation (20) 
in the sense of distribution. In fact, we notice that from the Lipschitz character of G, we have 

JJ y{G{ey) - G(0))g,dxdy < Ge JJ y'^q^dxdy -^0, as e 0, 

thanks to the estimate on the tail above. Finally, (20) implies that yG{0)qo + DtMqQ is constant a.e., 
and this constant should be 0 since from the estimate on the tail above, we have that yqo G L^. We 
conclude that qo vanishes except when yG{0) + DtM{S) = 0. By conservation of the mass, we deduce 
the expression (21) for po- 















5 Comments on physical background 

The form of the equation Q corresponds, for E. coli chemotaxis, to the formalism in the physical 
literature. We have simplified the notations for mathematical clarity and we explain now how to relate 
our notations to known biophysical quantities. Here we have used the same biological parameters as 

in dang]. 


• The methylation level M{x,t) at equilibrium is related to the extra-cellular attractant profile S, 
by a logarithmic dependency 

M = M{S) = m„+E^, with /„(S)=h,(i±^). 

The constants mo, Kj, Ka represent the basic methylation level, and the dissociation constants 
for inactive, respectively active, receptors. Numerical values are given by a = 1.7, mo = 1, 
Ki = 18.2/rM, Ka = 3mM. 

• The receptor activity a(m, S) depends on the intracellular methylation level m and the extra¬ 
cellular chemoattractant concentration S such that 

a = (l -|- exp(A^iii)) with E = —a{m — mo) -|- fo{S) = —a{m — M{S)). (27) 

The coefficient N = 6 represents the number of tightly coupled receptors. 


The intracellular dynamics and tumbling frequency are given by 


f{m-M{S))=F{a) = kR{l-a/ao), X{m - M{S)) = Z{a) = zq + t 

Oo 


where a(m, S) is the receptor activity defined in (27). The parameter kn is the methylation rate, 


ao is the receptor preferred activity which is such that /(ao) = 0 and f'{ao) < 0. For the tumbling 
frequency, zq, El, r represent the rotational diffusion, the Hill coefficient of flagellar motors 
response curve and the average run time respectively. All these parameters can be measured 


biologically, their values are given by kR = 0.01s ^ 
r = 0.8s, H = 10. 


0.0005s 


-1 


uq = 0.5, zo = 0.14s 


-1 


• Two kinds of noise can be observed in the signaling pathway for E. coli, one is from the external 
fluctuation of the ligand concentration and the other is the internal noise from the random 
receptor-methylation and demethylation event [6j. For small complexes, the effect on the activity 
from the external noise is negligible compared to the internal noise. 


We refer the readers to [laiiiiiiiiiH] and the references therein for the detailed physical meanings of 
these parameters. 

In E. coli chemotaxis, the adaptation time is faster than the system time, i.e. the rescaled constant 
kR is large. Let kR = j, where e is the ratio between the system time and the adaptation time. For 
example, when the system time is 1000s and the adaptation time is 100s, e = 0.1. 

Moreover, the Hill coefficient H is large which indicates that A(m — M{S)) varies fast with respect 
to m — M{S). Therefore, the scaling introduced in ([^ is satisfied by E. coli chemotaxis. We can use 

/(r) = 1—^(l-|-exp(—iVar)) 
ao 
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Therefore, the function G(-) used in pTj ) is given by 

G{r) = = —-(l —^(l + exp(—A^ar)) with G(0) = —/'(O) 

r r \ an ' 


Besides, from Theorem 2.1 for the case without noise, using y = \ yields 


T{DtM,v,v') = A{y,v,v') = zo+ T 


(1 + exp{-Na€y)) 


-H 


-.H 


= Zq + T 


-1 


1 + exp iyNaeDt 


-H 


a, 


And from Theorem |2.1[ when the noise in the methylation level is considered, 

V 2 


T{DtM,v,v') = 


GjO) 

2tt 


A{y,v,v')e 


G(0) I DtMy 

2 \y+ G(0) ) 


dy. 


Na 
4ao ' 


Since the run duration last longer when bacteria encounter an increasing gradient of chemotactic 
molecules, this lead to higher bacteria density at the place where the ligand concentration S is higher. 
This phenomena is well explained by the classical Keller-Segel model which can be considered as the 
parabolic limit of the kinetic-transport models. However, recent experimental observation shows that 
higher ligand concentration leading to higher bacteria density is only valid in a spatial-temporal slow- 
varying environment. When the ligand concentration varies fast, there exists a phase shift between 
the mass center of ligand concentration and of the bacterias [23j . This is due to the memory effect in 
the slow methylation adaptation rate. In the limiting kinetic model, DtM{S) takes into account the 
memory effect along the trajectory of the moving cell [20]. Then, two interesting questions come: are 
these two kinds of memory effect the same? Can we see the phase shift between the mass center of 
ligand concentration and of the bacteria in the limiting kinetic model? 


6 Numerical illustrations 

We performed numerical simulations using the method SPECS m- It is a cell based model that takes 
into account the evolution of each cell intracellular methylation level, which determines the tumbling 
frequency of each bacteria. As explained in [HldH], SPECS and the kinetic model that incorporates 
intracellular chemo-sensory system Q show a quantitative match. As in jl8l 119] . we choose the two 
velocity kinetic model and a periodic ID traveling wave concentration 

S{x,t) = S'o + S'a sin[^(x - ut)], 

which is spatial-temporal varying and where the wavelength i is equal to the length of the domain. 
We compare the numerical results of SPECS and the limiting kinetic model in Eigure Upwind 
scheme is used for the transport terms and periodic boundary conditions are considered. The density 
is scaled at the order of 10“^, it is the ratio between the actual cell number and the total number. In 
the Eigure]^ the x axis represents the remainder of x — ut mod i, i.e. we keep tracking the wave front 
of the periodic traveling wave. 

Two different wave velocities u = 0.4^m/s and u = 8ym/s are considered. We compare the density 
profiles p = fypdv and the cell flux J = JyVpdv. When the concentrated wave moves slowly, the 
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Figure 1: Numerical comparison between the limiting kinetic model and SPECS. The steady state 
profiles of density p (top) and current J (bottom) are presented. Left: u = 0.4^m/s; Right: u = 
8pm/s. In the subfigures, histograms and solid lines are from SPECS, while the dash star and dash 
dotted lines are calculated using the limiting kinetic equation. Parameters used here are [L]o = 500^M, 
[L]a = lOOpM, £ = 800pm. We have used 20, 000 cells for simulation with SPECS. 


limiting kinetic model gives good consistency, however in the fast-varying environment, the density 
and cell flux profiles are different for SPECS and the limiting kinetic model. We refer the readers to 
uadi for more detailed discussions and physical explanations. 

This numerical experiment shows that the memory effect using the model based on DtM(S) is 
different from the memory effect for the complete model (|^ when fast external chemoattractant waves 
are considered. This phenomena can be explained by the slow adaptation rate in the methylation 
level and the memory along the trajectory compared to the phase shift. In this fast wave regime, our 
mathematical results do not apply because the scaling assumptions are not satisfied. 
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